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Abstract 

Unknown quantum electronic conductance across nanojunctions made of silicon-doped carbon wires between 
carbon leads is investigated. This is done by an appropriate generalization of the phase field matching theory for the 
multi-scattering processes of electronic excitations at the nanojunction and the use of the tight-binding method. Our 
calculations of the electronic band structures for carbon, silicon, and diatomic silicon carbide are matched with the 
available corresponding density functional theory results to optimize the required tight-binding parameters. Silicon 
and carbon atoms are treated on the same footing by characterizing each with their corresponding orbitals. Several 
types of nanojunctions are analyzed to sample their behavior under different atomic configurations. We calculate for 
each nanojunction the individual contributions to the quantum conductance for the propagating a, tt, and a* 
electron incidents from the carbon leads. The calculated results show a number of remarkable features, which include 
the influence of the ordered periodic configurations of silicon-carbon pairs and the suppression of quantum 
conductance due to minimum substitutional disorder and artificially organized symmetry on these nanojunctions. 
Our results also demonstrate that the phase field matching theory is an efficient tool to treat the quantum 
conductance of complex molecular nanojunctions. 
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Background 

Quantitative analysis of electronic quantum transport in 
nanostructures is essential for the development of nano- 
electronic devices [1]. The monatomic linear carbon wire 
(MLCW) systems are expected in this context to have 
potentially interesting technological applications, in par- 
ticular as connecting junction elements between larger 
device components [2]. In this respect, electronic quan- 
tum transport properties are the key features of such wire 
nanojunctions [3]. 

Carbon exists in nature under a wide range of allotropic 
forms as the two-dimensional graphene [4], the cage 
fullerenes [5], and the quasi one-dimensional carbon nan- 
otubes [6] . These forms exhibit exceptional physical prop- 
erties and can be considered as promising components for 
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future nanodevices [7]. The discovery of MLCW, [8-14] 
turns the attention to another intriguing carbon allotropic 
form. In the experiment conducted recently by Jin et al 
[14], MLCW was produced by directly removing carbon 
atoms row by row from the graphene sheets, leading to a 
relatively stable freestanding nanostructure. 

At present, the available experimental data do not pro- 
vide essential knowledge about the electronic properties 
of MLCW systems, and only theoretical studies shed 
some light on these properties. Furthermore, although 
the MLCW systems were investigated for a long time 
from the theoretical point of view [15-26], their interest 
was not highlighted until recently due to the open atten- 
tion paid to other carbon allotropic forms. It has been 
shown in particular that from the structural point of view, 
MLCW can form either as cumulene wires (interatomic 
double bonds) or polyyne wires (alternating interatomic 
single and triple bonds) [14,17,19,27,28]. However, there 
is no straightforward answer as to which of these two 
structures is the favorable one; experimental studies do 
not give a satisfactory answer, and theoretical calculations 
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yield provisions which depend on applied computational 
methods. Density functional theory (DFT) calculations 
predict double-bond structures [29,30], whereas ab initio 
Hartree-Fock (HP) results favor alternating bond systems 
[15-18,27]. This situation arises from the fact that DFT 
tends to underestimate bond alternation (second-order 
Jahn-Teller effect), while HF overestimates it [27]. 

More recently, first-principle calculations have indi- 
cated [31] that both structures are stable and present 
mechanical characteristics of a purely one-dimensional 
nanomaterial. Moreover, on the basis of the first-principle 
calculations [31-42], the cumulene MLCW wires are 
expected to be almost perfect conductors, even bet- 
ter than linear gold wires [29], while the corresponding 
polyyne wires are semiconducting [41]. It is also worth 
noting that the MLCW cumulene system may exhibit con- 
ductance oscillations with the even and odd numbers of 
the wire atoms [28,42]. 

In the present work, we consider in particular the prob- 
lem of the electronic quantum transport across molecular 
nanojunctions made up of silicon-doped carbon wires, 
prepared in ordered or substitutionally disordered config- 
urations as in the schematic representation of Figure 1, 
where the nanojunctions are between pure MLCW wire 
leads. This problem has not been considered previously 
and is still unsolved to our knowledge. The interest in the 
quantum transport of such nanojunctions arises from the 
fact that chemical defects or substitutional disorder may 
have a significant impact on their transport properties 
[43]. Chemical impurities doping the nanojunction may 
even allow the control of the transport for such nanostruc- 
tures [44]. The properties of the nanoelectronic device 
and its functionality may hence be greatly affected or even 
built on such ordered and disordered configurations. The 



interest in silicon carbide, furthermore, stems from the 
fact that it is considered a good substrate material for 
the growth of graphene [45] and may produce interesting 
effects in its interactions with Si or C [46] . 

The electrons which contribute to transport present 
characteristic wavelengths comparable to the size of 
molecular nanojunctions, leading to quantum coherent 
effects. The transport properties of a given nanojunc- 
tion are then described in terms of the Landauer-Biittiker 
theory [47,48], which relates transmission scattering to 
quantum conductance. Several approaches have been 
developed in order to calculate the scattering trans- 
mission and reflection cross sections in nanostructures, 
where the most popular are based on first-principle cal- 
culations [49,50] and semiempirical methods using the 
non-equilibrium Greens function formalism [51,52]. 

In the present work, we investigate the electronic scat- 
tering processes on the basis of phase field matching 
theory (PFMT) [53,54], originally developed for the scat- 
tering of phonons and magnons in nanostructures [55-59]. 
Our theoretical method is based on appropriate phase 
matching of the Bloch states of ideal leads to the local 
states in the scattering region. In this approach, the 
electronic properties of the system are described in the 
framework of the tight-binding formalism (TB) which 
is widely exploited for electronic transport calculations 
[54,60-63] and for simulating the STM images of nanos- 
tructures [64,65]. In particular, we employ the appropriate 
Slater-Koster [66] type Hamiltonian parameters calcu- 
lated on the basis of the Harrisons tight-binding theory 
(HTBT) [67]. The PFMT method, which is formally equiv- 
alent to the method of non-equilibrium Green s functions 
[68], can be considered consequently as a transparent 
and efficient mathematical tool for the calculation of the 
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Figure 1 Schematic representation of finite silicon-doped carbon wire nanojunction between two semi-infinite quasi one-dimensional 
carbon leads. The irreducible region and matching domains are distinguished (please see subsection 'Phase field matching theory' in the 'Methods' 
section for more details). The binding energies for a given atomic site and the coupling terms between neighbor atoms with corresponding 
interatomic distances are depicted. The n and n' indices for the coupling parameters are dropped for simplicity. 
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electronic quantum transport properties for a wide range 
of molecular- sized nanoj unction systems. 

The present paper is organized in the following man- 
ner. In the 'Methods' section, we give the detailed dis- 
cussion of theoretical PFMT formalism. Our numerical 
results, which incorporate propagating and evanescent 
electronic states, are presented per individual lead modes 
in the 'Results and discussion' section. Also presented are 
the total conductance spectra; they are compared with 
results based on first-principle calculations when avail- 
able. Finally, the discussion and conclusions are given in 
the 'Conclusions' section. Appropriate appendices which 
supplement the theoretical model are also presented. 

Methods 

Theoretical model and propagating states 

The schematic representation of the system under study 
with an arbitrary nanoj unction region is presented in 
Figure 1. With reference to the Landauer-Biittiker the- 
ory for the analysis of the electronic scattering processes 
[47,48], this system is divided into three main parts, 
namely the finite silicon-doped carbon wire nanojunc- 
tion region, made up of a given composition of carbon 
(black) and silicon (orange) atoms, and two other regions 
to the right and left of the nanoj unction which are semi- 
infinite quasi one-dimensional carbon leads. Moreover, 
for the purpose of quantum conductance calculations, the 
so-called irreducible region and the matching domains 
are depicted (see the 'Phase field matching theory' sub- 
section for more details). Figure 1 is used throughout the 
'Methods' section as a graphical reference for analytical 
discussion. 

The system presented in Figure 1 is described by the 
general tight-binding Hamiltonian block matrix: 



H = 



H 



0 
0 



(1) 



0 Ha^,A^_i Em,N ^''N+hN 

0 0 Hiv+i,Ar En-\-i,n-\-i 
0 0 



This is defined in general for a system of Nx inequiv- 
alent atoms per unit cell, where Ni denotes the number 
of basis orbitals per atomic site, assuming spin degen- 
eracy. In Equation 1, Eij denotes on-diagonal matrices 

composed of both diagonal s'^'^ and off-diagonal h^'l^^ 
elements for a selected unit cell. In contrast, the H/,y matri- 
ces contain only off-diagonal elements for interactions 
between different unit cells. The index a identifies the 
atom type, C or Si, on the nth. site in a unit cell. Each 
diagonal element is characterized by the lower index / for 



the angular momentum state. The off-diagonal elements 
describe the m-type bond, {m = a, n), between / 



and V nearest-neighbor states. The index P identifies the 
types of interacting neighbors, C-C, Si-Si, or Si-C. 



The h 



elements are consistent with the Slater- 



1,1', m 

Koster convention [66] and may be expressed in the 
framework of the HTBT [67] by the following: 



h 



l,l',m 



11,1' 



(2) 



where values are the dimensionless Harrison coef- 

ficients; Me, the electron mass in vacuum; and dp, the 
interatomic distance for interacting neighbors. Explicit 
forms of the E/,y and Hij matrices are given in Appendix 
1. The tight-binding parameter schemes are illustrated 
in Figure 1; however, it is noteworthy that the n and 
indices for coupling parameters are dropped for simplicity 
in this figure. 

In our calculations, the single-particle electronic wave 
functions are expanded in the orthonormal basis of local 
atomic wave functions 0/(r) as follows: 



vl/(r,k) = ^c/(r« 

l,n,N 



RiV,k)0/(r-Riv,k). 



(3) 



In Equation 3, k is the real wave vector; R^v, the position 
vector of the selected unit cell; and r^, the position vector 
of the nt\\ atom in the selected unit cell. For ideal leads, the 
wave function coefficients c:/(r^ — Ra/^, k) are characterized 
under the Bloch-Floquet theorem in consecutive unit cells 
by the following phase relation: 



CiiVn - RiV+1, k) = zciixn - RiV, k), 
where z is the phase factor 



(4) 



(5) 



which corresponds here to waves propagating to the right 
(+) or to the left (-). 

The electronic equations of motion for a leads unit cell, 
independent of N, may be expressed in a square matrix 
form, with an orthonormal minimal basis set of local wave 
functions as follows: 



(El-Md) xc(k,£) = 0. 



(6) 



E stands for the electron eigenvalues, and I is the iden- 
tity matrix, while the dynamical matrix contains the 
Hamiltonian matrix elements and the z phase factors; 
c(kfE) is theNx x Ni size vector defined as follows: 



c(k,£) 



C5(ri,k,£) 
c^^(ri,k,£) 



'c/(ri,k,£) ■ 
^c/(r„,k,£)_ 



(7) 
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Equation 6 gives the Nx x A// eigenvalues with cor- 
responding eigenvectors which determine the electronic 
structure of the lead system, where / under the vector c/ 
corresponds to Ni = 4 orbitals s.px.Py.Pz- Note that the 
choice of an orthonormal minimal basis set of local wave- 
functions may result initially in an inadequate description 
of the considered electronic eigenvalues. However, as can 
be seen later, the proper choice of the TB on-site ener- 
gies and coupling terms allows us to to obtain agreement 
with the DFT results. This is a systematic procedure in our 
calculations. 

Evanescent states 

The complete description of electronic states on the ideal 
leads requires a full understanding of the propagating 
and evanescent electronic states on the leads. This arises 
because the silicon-doped nanojunction breaks the per- 
fect periodicity of the infinite leads and forbids a formula- 
tion of the problem only in terms of the pure Bloch states 
as given in Equation 5. Depending on the complexity of a 
given electronic state, it follows that the evanescent waves 
may be defined by the phase factors for a purely imaginary 
wave vectors k = m such that 



z = z±= e^'^''", 



or for complex wave vectors k = ki -\- iK2 such that 

Z = Z±= 6^^'"'-"^^'^ 



(8) 



(9) 



The phase factors of Equations 8 and 9 correspond to 
pairs of hermitian evanescent and divergent solutions on 
the leads. Only the evanescent states are physically con- 
sidered where spatial evanescence occurs to the right and 
left, away from the nanojunction localized states. It is 
important to note that the /-type evanescent state corre- 
sponds to energies beyond the propagating band structure 
for this state. 

The functional behavior of z(E) for the propagating and 
evanescent states on the leads may be obtained by various 
techniques. An elegant method presented previously for 
phonon and magnon excitations [59] is adapted here for 
the electrons. It is described on the basis of Equations 4 
and 6 by the generalized eigenvalue problem for z: 



El 



- Em,n ^n,n- 
I 0 



— z 



-H 



t 

■N,N- 

0 



0 



c(Ra/',z,£) 



0. 



(10) 



Equation 10 gives the 2NxNi eigenvalues as an ensem- 
ble ofNxNi pairs of ;^ and z~^. Only solutions with \z\ = 1 
(propagating waves) and \z\ < 1 (evanescent waves) 
are retained as physical ones. In Equation 10, k is then 
replaced by the appropriate energy E variable. Further- 
more, for systems with more than one atom per unit cell. 



the matrices H^^m-i and h]^^_^ in this procedure are 
singular. In order to obtain the physical solutions, the 
eigenvalue problem of Equation 10 is reduced from the 
2NxNi size problem to the appropriate 2Ni one, using the 
partitioning technique (please see Appendix 2). 

Phase field matching theory 

The scattering problem at the nanojunction is considered 
next. An electron incident along the leads has a given 
energy E and wave vector k, where E = Ey (k) denotes the 
available dispersion curves for 7 = 1, 2,..., F propagating 
eigenmodes, where F corresponds to the total number of 
allowed solutions for the eigenvalue problem of phase fac- 
tors in Equation 10. In any given energy interval, however, 
these may be evanescent or propagating eigenmodes and 
together constitute a complete set of available channels 
necessary for the scattering analysis. 

The irreducible domain of atomic sites for the scatter- 
ing problem includes the nanojunction domain itself, (N G 
[ 0, D — 1]), and the atomic sites on the left and right leads 
which interact with the nanojunction, as in Figure 1. This 
constitutes a necessary and sufficient region for our con- 
siderations, i.e., any supplementary atoms from the leads 
included in the calculations do not change the final results. 
The scattering at the boundary yields then the coherent 
reflected and transmitted fields, and in order to calculate 
these, we establish the system of equations of motion for 
the atomic sites {N e[—hD])of the irreducible nanojunc- 
tion domain. 

This procedure leads to the following general matrix 
equation: 



0. 



(11) 



^nano is a (D + 2) X (Z) + 4) matrix composed of the block 
matrices (El — Em,n — ^n,n-i — ^mm-i^' ^^^^^ 
vector V of dimension £) + 4 is given as follows: 

c/(ri -R_2,£) 



V = 



c/(r« -R_2,£) 



c/(ri -Ri)+i,£) 



(12) 



_c/(r„ -Ri)+i,r) 

Since the number of unknown coefficients in 
Equation 11 is always greater than the number of 
equations, such a set of equations cannot be solved 
directly. 

Assuming that the incoming electron wave propagates 
from left to right in the eigenmode y over the interval 
of energies E = Ey^ the field coefficients on the left and 
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right sides of the irreducible nanojunction domain may be 
written as follows: 

~ Zy, Ey) = CliVyi, Zy , Ey)Zy^ 

r 

+ ^ci(rn,Zyf,Ey)z^fryy(Ey) 



The rectangular sparse matrix in Equation 15 has the 
(£) + 4) X (D -\- 2) size. The vectors r^^^/ and ty^yf are 
column vectors of the backscattering and transmission 
Hilbert basis. 

Substituting Equation 15 into Equation 11 yields an 
inhomogeneous system of equations as follows: 



for ^ -1, 



(13) 



Cf - RNyZy,Ey) = y^^Cl(rnyZy'yEy)Z^fty^y^(Ey) 

y' 

for N^D, (14) 

where e F is an arbitrary channel into which the inci- 
dent electron wave scatters, and ciixn, Zy.Ey) denotes the 
the eigenvector of the lead dynamical matrix of Equation 6 
for the inequivalent site n at Zy and Ey . The terms Vy^yf 
and ty^y' denote the scattering amplitudes for backscat- 
tering and transmission, respectively, from the y into the 
y^ eigenmodes and constitute the basis of the Hilbert 
space which describes the reflection and transmission 
processes. 

Equations 13 and 14 are next used to transform the 
(D + 2) X (D + 4) matrix of the system of equations of 
motion. Equation 11, into an inhomogeneous (D + 2) x 
iP + 2) matrix for the scattering problem. This procedure 
leads to the new form of the following vector: 



z 0 
0 1 



1 0 

0 z 

0 z^ . 



c/(ri -Ro,£j/) 



c/(ri -R£)_i,£j,) 
c/(r„ - RD_i,£y) 

ty,y' 



Cl(ri,Zy,Ey)Zy^ 



Cl(rn,Zy,Ey)Zy^ 

ci(ri,Zy,Ey)z-^ 



Cl{Yn,Zy,Ey)Zy^ 

0 



Mx 



c/(ri 



-Ro,^k) 



(16) 



c/(r« -Ro,£y) 



c/(ri -Ri)_i,£y) 



c/(r„ -Ro-i^Ey) 

ty^y' 



In Equation 16, M is the matched (D + 2) x (D + 2) 
square matrix, and the vector of dimension {D + 2) which 
incorporates the M^" and elements, regroups the 
inhomogeneous terms of the incident wave. The explicit 
forms of the M matrix elements and and MJ^ vectors are 
presented in Appendix 3. 

In practice, Equation 16 can be solved using standard 
numerical procedures, over the entire range of available 
electronic energies, yielding the coefficient c/ for atomic 
sites on the nanojunction domain itself as well as the T 
reflection ry^yf{E) and the T transmission ty^yf{E) coeffi- 
cients. 

The reflection and transmission coefficients give the 
reflection Ry^y'iE) and transmission Ty^y'iE) probabili- 
ties, respectively, by normalizing with respect to their 
group velocities Vy in order to obtain the unitarity of the 
scattering matrix as follows: 



RyyiE) 



'y,y 



-y,y 



(E) = ^\tyy(E)\\ 



(17) 



(18) 



(15) 



where Vy = Vy(E) denotes the group velocity of the 
incident electron wave in the eigenmode y. The group 
velocities are calculated by a straightforward procedure 
as in Appendix 4. For evanescent eigenmodes, v^/ = 0. 
Although the evanescent eigenmodes do not contribute to 
the electronic transport, they are required for the com- 
plete description of the scattering processes. 

Furthermore, using Equations 17 and 18, the over- 
all reflection probability, Ry{E), for an electron incident 
in the y eigenmode and the total electronic reflection 
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probability, R(E), from all the eigenmodes may be ex- 
pressed, respectively, as follows: 

r r 
Ry(E) = J2^yy(E) and R(E) = J2^y(E). (19) 

y 

Similarly, for transmission probabilities, we may write 
the equivalent equations as follows: 

r r 
Ty(E) = Tyy(E) and T(E) = ^ Ty(E), (20) 

y 

The Ty (E) and T(E) probabilities are very important for 
the electronic scattering processes since they correspond 
directly to the experimentally measurable observables. 
Likewise, the total transmission T(Ey) allows to calcu- 
late the overall electronic conductance. In this work, we 
assume the zero-bias limit and write the total conductance 
in the following way: 

G(Ef) = GoT(Ef), (21) 

In Equation 21, Go is the conductance quantum and 
equals 2e^/h, Due to the Fermi-Dirac distribution, G{Ep) 
is calculated at the Fermi level of the perfect lead band 
structure since electrons only at this level give the impor- 
tant contribution to the electronic conductance. The 
Fermi energy can be determined using various methods 
where, in the present work, Ep is calculated as the basis of 
the density of state calculations. 

Results and discussion 

The tight-binding model and basic electronic properties 

In this section, we present the results of our model cal- 
culations for the electronic structure of carbon, silicon, 
and silicon carbide wires under study. Our results are vali- 
dated by comparison with DFT calculations [29,69], which 
allow us to establish unambiguously our choice of the 
tight-binding parameters for these systems. 

In principle, we can develop our model calculations for 
the nanojunctions and their leads using any adequate type 
of orbitals; even a single orbital suffices to calculate the 
electronic quantum transport for carbon nanojunctions 
[44]. However, this approximation is inadequate for silicon 
atoms. To treat both types of atoms on the same footing, 
we thus characterize the atoms by the electronic states 2s 
and 2p for carbon and by 35 and 3p for silicon. Such a 
scheme gives us four different orbitals, namely 5, py, 
and pz, for both types of atoms. 

In the present work, our TB parameters are effectively 
rescaled from the Harrisons data in order to match our 
model calculations for the electronic structure with those 
given by the DFT. The utilized TB parameters are pre- 
sented in Table 1 in comparison with the values given by 



Table 1 Tight-binding parameters and Harrison's 
dimensionless coefficients proposed in this work and 
compared with original values 

Harrison TB parameters Present TB parameters 

a 

£s 
£p 

hs,p,a 
^P,P,(y 

hp,p,7t 

± 

The values of the tight-binding parameters, e"'" and h^'i,'^ (in eV), and 
Harrison's dimensionless coefficients, rm',mi as proposed in this work and 
compared with the original values by Harrison [67]. Please note that the 
distance-dependent h^'i,f parameters are computed for the appropriate 
interatomic spacings (in A), tabulated below, and assumed after the works of 
Tongay et al. [29] and Bekaroglu et al. [69]. In order to keep table transparent 
indices, n and n' for e"'" and h^'^ are omitted. 



Harrison. It is worthy to note that the values of the on-site 
Hamiltonian matrix elements s^'^ are identical for states 

Px>Pyy and pz. The off-diagonal distance-dependent hj'l!^ 
elements are calculated on the basis of Equation 2. For 
symmetry considerations, these latter elements are posi- 
tive or negative, also hs,p,a = ^s,p,a = 0 and /z^,^,^ = 77^,^,^ = 
0, for py and pz, and /z^,^,;^ = hpy^py^jj: = hp^^p^^jj; = hp^^p^^jj; = 
0 [70] . Table 1 is supplemented for the reader by Figure 2 
which gives the dependence of the hopping integrals with 
distance as calculated in the present paper (continuous 
curves), in comparison with the Harrisons data (open 
symbols). 

Figure 2 clearly indicates the fact that qualitatively, both 
Harrison's and our rescaled coupling parameters for sil- 
icon, carbon and diatomic silicon carbide wires, present 
the same functional behavior, confirming the desired con- 
servation of their physical character. However, most of 
the rescaled coupling parameters have somehow smaller 
values than those initially proposed by Harrison; this 
trend can be also traced in Table 1 for the onsite param- 
eters. This difference stems from the influence of the 
low-coordinated systems are considered here, whereas 
the initial Harrison values are given to match tetrahedral 
phases [67]. Another general observation can be made for 
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Figure 2 The nearest-neighbor tight-binding coupling parameters with the interatomic distance (A, B, C). The curves represent our 
calculated TB results in comparison with those calculated using the Harrison parameters (squares, triangles, circles). 



tlie tiglit-binding parameters of the a -type interactions 
{the hs,p,cj and hp^p^^ ones), wliicli present vcmch closer val- 
ues over the considered interatomic distance range than 
in the case of Harrison s data. 

Our calculated electronic band structures for silicon, 
carbon, and diatomic silicon carbide infinite wires (con- 
tinuous curves) are presented in Figure 3 in compari- 
son with the DFT results [29,69] as in the right-hand 
side of the figures. We note for the carbon and silicon 
structures that our TB parameters correctly reproduce 
the DFT results up to energies slightly above the Fermi 
level. Electronic branches in the regions of high energies 
are in qualitative agreement. In the case of the diatomic 
silicon carbide structure, some of the electronic states 
perfectly match the DFT results even for high-energy 
domains. The left-hand side of Figure 3 compares our 
results (continuous curves) with those from the older 
TB values given by Harrison (open symbols); as is seen, 
our TB parameters constitute the most optimal set 
for the electronic transport calculations since their cor- 
responding electronic band structures conform to the 
appropriate energy ranges highlighted by the DFT results 
and, what is even more important, correctly reproduce the 
Fermi level. 

In Figure 3A,B for silicon and carbon, the red and blue 
colors correspond, respectively, to the a and a* bands. 
These arise from the spx orbital hybrids where the lowest 
lying bands are always occupied by two electrons. Bands 
marked by the red color have the n character and are 



degenerate. Their origin in the py and pz orbitals allows 
them to hold up to four electrons. In Figure 3C for the 
diatomic silicon carbide, starting from the band structure 
minimum, consecutive bands have their origin in the fol- 
lowing orbitals: carbon 35 (red band), silicon 35 (green 
band), carbon ?>p (blue and black bands), and silicon ?>p 
(orange and violet bands). The blue and orange colors 
for the silicon carbide electronic structure indicate two 
doubly degenerate tt -type bands. 

The metallic or insulating character of the considered 
atomic wires, following the Fermi level, is appropriate only 
when the wires are infinite. It is well known that this char- 
acter can change for the case of finite size wires with a 
limited number of atoms or due to the type and quality of 
the leads. 

Numerical characteristics for the carbon leads 

In general, the infinite carbon wires which are consid- 
ered as the leads in our work, present electronic band 
structure characteristics which incorporate not only prop- 
agating (see Figure 3A), but also evanescent states. Both 
of these types of states, which are derivable from the gen- 
eralized eigenvalue problem as presented in Equation 10, 
constitute a complete set over the allowed energies for the 
electrons incident along the leads, which can be further 
scattered at the considered nanojunction. This complete 
set of eigenstates is used as the basis for the numerical 
calculations of the quantum conductance presented in the 
'Transport properties' subsection. 
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Figure 3 Electronic structures of carbon (A), silicon (B), and diatomic silicon carbide (C). These structures are for infinite linear atomic wires 
presented over tine first Brillouin zone (p = kd g[ -tt, tt]. Our calculated results (continuous curves), represented by a color scheme (details in the 
text), are compared on the right-hand side with the first-principle results (closed circles, e[0,7i]) [29,69] and on the left-hand side with results 
calculated using Harrison TB parameters [67] (diamonds, (p g[ -tt, 0]). Our calculated Fermi levels are given as the zero-reference energies, and the 
calculated electronic DOS in arbitrary units are presented in the right-hand column. 



Figure 4A presents the three-dimensional representa- 
tion of the solutions of Equation 10 as a set of generalized 
functionals z(E) for the <7, cr*, and n electronic states of 
the carbon leads. As described by Equations 5, 8, and 9, 
the eigenstates in Figure 4A characterized by |z| = 1 cor- 
respond to the propagating electronic waves described by 
the real wave vectors, whereas those by |:^| < 1 corre- 
spond to the evanescent and divergent eigenstates for the 
complex wave vectors. Furthermore, for convenience, the 
corresponding moduli of the complex z factors are pre- 
sented in Figure 4B. Note that \z\ = 1 solutions may 
be grouped into pairs for the two directions of propaga- 
tion linked by time-reversal symmetry. Due to the fact 



that each of these two solutions provides the same infor- 
mation, we consider waves propagating only from left to 
right. However, this is not true for the \z\ < 1 solutions 
which are always considered for both left and right as 
spatially evanescent. As can be seen in Figure 4, the gen- 
eralized results for a, a*, and tt states are represented by 
the same colors as the corresponding states in Figure 3A, 
following their propagating character for |z| = 1, and 
further extended to the physically \z\ < 1 evanescent 
solutions. 

Figure 4 provides a more complete description for the 
electronic states of a given system compared to a typical 
band structure representation as in Figure 3, since both 
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the propagating and evanescent states are shown. Such 
a general representation clearly indicates the importance 
of the evanescent eigenstates for a full description of the 
scattering problem presented in the 'Transport properties' 
subsection. The energies considered in our calculations 
correspond to the range within the band structure bound- 
aries, marked by two vertical dotted lines in Figure 4B. 
As a consequence, not only the propagating states, but 
also the evanescent solutions are included in the quantum 
conductance calculations in the 'Transport properties' 
subsection. 

Transport properties 

In this subsection, the electronic transport properties 
of nanojunction systems composed of silicon-doped car- 
bon wires between carbon leads are calculated using the 
PFMT method. Figure 5A presents a number of these sys- 
tems where we indicate the irreducible domains by the 
shaded grey areas. Note that these systems are always 
composed of finite nanojunction regions of silicon and 
carbon atoms, coupled with two carbon semi-infinite 
leads. The first three systems of Figure 5 correspond 
to periodic diatomic silicon carbide nanojunctions com- 
posed of 1, 2, and 3 Si-C atomic pairs, respectively. The 
next system corresponds to a nanojunction with a sub- 
stitutional disorder, composed of three carbon and three 
silicon atoms. The last is a symmetric nanojunction of five 
silicon atoms and only one carbon atom in the middle. 
Figure 5B presents the group velocities of electrons in the 
carbon leads. 

The calculated transmission and reflection scattering 
cross sections for each of the four available transport 
channels are presented in Figure 6. Each row of the 
figure corresponds to a nanojunction system (NS) as 



follows: Figure 6A,B,C for NS 1, Figure 6D,E,F for NS 
2, Figure 6G,H,I for NS 3, Figure 6J,K,L for NS 4, and 
Figure 6M,N,0 for NS 5. The red and green continuous 
curves represent the transmission and reflection spec- 
tra, respectively. The blue histograms correspond to the 
free electronic transport on the carbon leads, i.e., to the 
electronic transport on the perfect infinite quasi one- 
dimensional carbon wire over the different propagating 
states. These histograms constitute the reference to the 
unitarity condition which is used systematically as a check 
on the numerical results. The leads' Fermi level is marked 
by a dashed line and set as a zero-energy reference. Under 
the zero-bias limit, the total conductance is calculated at 
this Fermi level. 

In Figure 6, the transmission spectra present strong 
scattering resonances, showing an increasing complex- 
ity with the increasing size and configurational order of 
the nanojunctions. The valence a state exhibits negligible 
transmission for all of the considered nanojunctions. The 
degenerate tt states and the a* state present in contrast 
the finite transmission spectra. However, it is only the tt 
states which cross the Fermi level, giving rise to electronic 
conductance across the nanojunction within the zero-bias 
limit. 

In particular, the first three considered systems rep- 
resent increasing lengths of the diatomic silicon carbide 
nanojunction with the increasing number of ordered Si- 
C atomic pairs. The transmission at the Fermi level for 
these systems is nonzero (see Figure 3C), which contrasts 
with the insulating character of the infinite silicon car- 
bide wire. One can connect this finite transmission to 
the indirect bandgap (A) around the Fermi level for the 
diatomic silicon-carbide infinite wire (for more details, 
please see Figure 3C). This gap, A ~ 1.5 eV, is indeed 
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Figure 5 Schematic representation of the five nanojunction 
systems and group velocities for propagating band structure 
modes. (A) Schematic representation of tine five nanojunction 
systems composed of silicon and carbon atoms between 
one-dimensional carbon leads considered in the present work. The 
irreducible domains are marked by the shaded grey areas, whereas 
for the other cases, only the irreducible domains are shown. (B) The 
group velocities for the propagating band structure modes on the 
carbon leads. 



related to the difference between the binding energies 
of the silicon and carbon atoms and corresponds to an 
effective potential barrier for the propagating tt- state elec- 
trons. As the wire length increases by adding Si-C atomic 
pairs, as for systems 1 to 3 of Figure 5B, the transmis- 
sion decreases due to cumulative barrier effects. We note 



that a similar effect for the monovalent diatomic copper- 
cobalt wire nanojunctions has been observed in a previous 
work [54]. 

Furthermore, it is instructive to compare the scatter- 
ing spectra for the degenerate tt states, for nanojunction 
systems 3 and 4. These two systems contain identical 
numbers of silicon and carbon atoms; however, system 
3 is an ordered configuration of Si-C pairs, whereas sys- 
tem 4 presents substitutional disorder of the atoms. It is 
seen that the disorder suppresses the conductance of the 
TT -state electrons at the Fermi level within the zero-bias 
limit. Another general observation can be made from the 
results for nanojunction system 5 which contains more 
silicon than carbon atoms. Despite the finite size of this 
system, which is comparable to system 4, and despite 
the structural symmetry of its atomic configuration, the 
electronic transmission is suppressed at the Fermi level 
within the zero-bias limit. This implies that one of the 
main observations of our paper is that structural sym- 
metry on the nanojunction is not a guarantee for finite 
transmission in the case of the multivalence diatomic wire 
nanojunctions. 

Figure 6 also shows that the transmission spectra for 
the a* state are close to unity over a significant range of 
energies from approximately 1 to 7 eV for all of the five 
nanojunction systems. This result may prove useful for 
the electronic conductance across silicon-doped carbon 
nanojunctions under finite bias voltages. 

In Figure 7, we present the total electronic conductance 
G(E) as a function of energy E and in units of Go = 2e^/h. 
for the considered nanojunction systems of a given length 
as depicted in Figure 5 (red). Moreover, the perfect elec- 
tronic conductance on the carbon leads (blue) is given in 
comparison and constitutes effectively the conductance 
of the infinite and perfect quasi one-dimensional carbon 
wire. In Figure 7, the Fermi level is indicated by the dashed 
line as a zero-reference energy, and G(E) is calculated 
from all the contributing eigenstates of Figure 6, including 
the two degenerate tt states. 

We note that the conclusions given for the results pre- 
sented in Figure 6 are also followed by the more gen- 
eral representation of the electronic transport depicted in 
Figure 7. Furthermore, the results presented in Figure 7 
confirm that only the electrons incident from the leads 
in the tt states are responsible for the electronic conduc- 
tance at the zero-bias limit, which is readable from the 
Fermi level position. However, for all considered systems, 
the conductance at the Fermi level is theoretically limited 
to the value of 2 Go, and the biggest conductance max- 
ima close to the perfect infinite carbon wire value of 3 Go 
can be observed only in the energy interval from approx- 
imately 1 to 7 eV hence for energies above the Fermi 
level. Once again, this follows our previous observations 
for the transmission results for the tt states concluded 
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Figure 6 Transmission and reflection probabilities across five types of silicon-doped carbon wires between two semi-infinite 
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from Figure 6. Nonetheless, only on the basis of the results 
presented in Figure 7 can we note that due to the summa- 
tion over all possible state contributions which constitute 
the G(E) spectra, not only the a* -state electrons, but also 
some of those in the degenerate n states contribute to the 
high conductance values in the cited energy intervals. This 
important observation proves that the a*- and tt -state 



electrons are of crucial importance for both the zero-bias 
quantum conductance of the silicon-doped carbon wires 
and the possible finite bias ones. This implies that the use 
of only a single orbital for the description of the carbon 
atoms will result in an inadequate description of the trans- 
port processes across low-coordinated systems containing 
these atoms. 
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Conclusions 

In the present work, the unknown properties of the quan- 
tum electronic conductance for nanojunctions made of 
siUcon-doped carbon wires between carbon leads are 
studied in depth. This is done using the phase field match- 
ing theory and the tight-binding method. The local basis 
for the electronic wave functions is assumed to be com- 
posed of four different atomic orbitals for silicon and 
carbon, namely the 5, px> Py> and pz states. 

In the first step, we calculate the electronic band 
structures for three nanomaterials, namely the one- 
dimensional infinite wires of silicon, carbon, and diatomic 



silicon carbide. This permits a matching comparison with 
the available corresponding DFT results, with the objec- 
tive to select the optimal TB parameters for the three 
nanomaterials. 

This optimal set of the tight-binding parameters is then 
used to calculate the electronic conductance across the 
silicon-doped carbon wire nanojunctions. Five different 
nanojunction cases are analyzed to sample their behav- 
ior under different atomic configurations. We show that 
despite the nonconducting character of the infinite silicon 
carbide wires, its finite implementation as nanojunctions 
exhibit a finite conductance. This outcome is explained 
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by the energy difference between the binding energies of 
the silicon and carbon atoms, which correspond to an 
effective potential barrier for the degenerate tt -state elec- 
trons transmitted across the nanojunction under zero- 
bias field. 

The conductance effects that may arise due to minimal 
substitutional disorder and to artificially organize symme- 
try considerations on the silicon carbide wire nanojunc- 
tion are also investigated. By exchanging the positions of 
two silicon and carbon atoms on an initial nanojunction to 
generate a substitutional disorder, we show that the total 
quantum conductance is suppressed at the Fermi level. 
This is in sharp contrast with the finite and significant 
conductance for the initial atomically ordered nanojunc- 
tion with periodic configurations of the silicon and carbon 
atoms. Also, the analysis of a silicon carbide nanojunc- 
tion of a comparable size as the one above, presenting 
symmetry properties, shows that quantum conductance is 
suppressed at the Fermi level. 

In summary, we note that the biggest maxima of 
the conductance spectra for the zero-bias limit can be 
observed for high energies for all of the considered sys- 
tems. This conclusion reveals the fact that electrons inci- 
dent from the leads in both a* and tt states are crucial for 
the considerations of the electronic transport properties 
of the silicon-doped carbon wire nanoj unctions. 



Appendix 1 

Explicit forms of the E' j and H' j matrices 

The explicit forms of the submatrices of Equation 1 are 
given in the following manner: 



where 



^1 ^2,1 
h2,l €2 



1 



t 

n,n—l 



(22) 



and 



^/ J 



(24) 



and 



'^s,l',m 



-l,m 



h 



l-l,l',m 



(25) 



Equations 22 and 23 denote NxNi square matrices, 
where matrix (Equation 23) is upper triangular. In this 
manner, component matrices (Equations 24 and 25) are of 
the dimension Ni x Ni, Additionally, matrix Si^jf always 
denotes diagonal matrix, while h//^y/ matrix is much more 
complex, with possible nonzero elements at every posi- 
tion. Please note that some of the h^'n elements can van- 
ish due to symmetry conditions and simplify the notation 
of the hifjf matrix. 



Appendix 2 

Partitioning technique 

The partitioning technique is a suitable method which 
allows to avoid the singularity problem of the H7v,a/^-i and 
^I/^A/^-i i^atrices and calculates only nontrivial solutions 
of Equation 10. Detailed discussion of the partitioning 
technique is presented in the work of Khomyakov and 
Brocks [71], and this section gives only our short remarks 
on this method. 

Following studies from Khomyakov and Brocks [71], 
Equation 10 is partitioned into two parts of Di — D2 and 
D2 sizes where 



0 • • • hi 



0 



0 0 
0 0 



• • hi,„_i hi,^ 



(23) 



and 



D2=NnNi. 



(26) 



(27) 



In Equation 27, parameter Nn stands for the order of 
nearest-neighbor interactions assumed in calculations, 
e.g., Nn = 1 for the first nearest-neighbor interactions. 
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On the basis of Equations 26 and 27, the reduced 2Ni 
eigenvalue problem is written as follows: 



LC2( 







Bi,i Bi,2 




0 


— z 


0 l2,2_ 





= 0. 



(28) 



At this point, we correct the misprint from the study of 
Khomyakov and Brocks [71] and write the submatrices of 
Equation 28 in the following form: 

Ai,i = £l2,2 - E2,2 - E2,i [£Ii,i - Ei,i]"^ Ei,2, (29) 
Ai,2 = -H2,2 - E2,i [£Ii,i - Ei,i]"' Hi,2, (30) 
Bi,i = 2 + h{ 2 ~ Ei,i] ^ Ei,2, (31) 

Bi,2 = Hi,2 - Ei,i]"' Hi,2. (32) 

Please note that the reduced problem of Equation 28 
gives 2Ni eigenvalues with 2Ni corresponding eigenvec- 
tors; this Nx times less than can be expected from a 
physical point of view. Nevertheless, those solutions can 
be easily separated into A/^A// eigenvalues and A/^A// eigen- 
vectors of a purely physical character. 

Appendix 3 

Explicit forms of the Mg, M!", and M2 components 

The submatrices of the matched (D + 2) x (D + 2) square 
matrix M in Equation 16, for a given / and ; indices, are 
given as follows: 



M/,y =El- Ei^i for 



Mij = -Hi^i-i for 



I =J 
D> i> I 
D>J>1 

i>2 , 
j = i-l 



(33) 



(34) 



M^^ = -H 



for ^ / < D + 1 , (35) 
j = i + l 

except for the submatrices which describe the boundary 
atoms of the system and those that are expressed in the 
following manner: 

Mi,i = -H-i-2Ci{Tn,ZyyEy)z^y 

+ (£1 — E_l,_l) Cl{Tn,Zy'-,Ey)Zy'-, 
M2,l = —\iQ-lCl{Vn,Zy',Ey)Zy', 

= -^r>-hD^l^^nyZyfyEy)z^,, 



M/)4.2,D+2 = -liD,D-\-lCl(rn,Zyf,Ey)z^^^ 

+ (£1 - Ed,d) Cl{Tn.Zy',Ey)^t. 



Finally, the M^^ and 'bA^ of Equation 16 vector compo- 
nents are written as follows: 



Mf = \l-i,-2Cl(Xn.Zy.Ey)zf 

+ (£1 - E_i,_i) ci{vn.Zy,Ey)Zy^, (37) 
M^^ = -Ho,-iCi(rn. Zy,Ey)z-\ (38) 



and 



Appendix 4 

Group velocities 

As specified in the 'Phase field matching theory sub- 
section, the group velocities for individual states can be 
calculated on the basis of Equation 6 rewritten in the 
following manner: 



[vI-V]y(k,£) = 0, 



(39) 



where v denotes the eigenvalues of Equation 39 which 
yields all required electron group velocities for each prop- 
agating state. Further, V is the Nx x Ni size matrix of the 
following form: 



V = 



dk ' 



(40) 



(36) 



Finally, v(Ra/^, k) stands for eigenvectors of the problem 
of Equation 39. We note that, usually, Equation 40 includes 
the constant part <i^/h, where h is the Planck constant. 
However, for the purpose of electronic conductance cal- 
culations within the PFMT approach, this term can be 
omitted due to the fact that only the ratios of the given 
group velocities are important (please see Equations 17 
and 18). 
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